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We report experimental evidence of electrostatic inhibition of fast electrons, generated in a highly 
resistive material upon irradiation with an intense ultra-short (10^^ W/cm^ , 100 fmsec) laser pulse. 
The experiment involves measurement of temporal evolution of self-generated magnetic pulses using 
pump-probe polarimetry. A comparison is made between the temporal behaviour of magnetic pulses 
generated with Aluminum and Glass targets. It is found that in contrast to Aluminium, self- 
generated magnetic pulse decays much faster in glass. This is attributed to the absence of return 
shielding currents in glass, which results in build up of electrostatic field, which in turn inhibits the 
movement of fast electrons. Fitting of experimental measurements using a one dimensional model, 
yields estimate of conductivity of Aluminium and glass, and penetration depth of hot electrons in 
these materials. 

PACS numbers: 52.38.Fz, 52.70.Ds, 52.70.Kz, 52.65.Rr 
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Interaction of intense ultrashort laser pulses with a 
solid leads to generation of large magnetic fields. Mag- 
netic fields up to gigagauss magnitudes have been pre- 
dicted in the overdense region of solid density plasmas 
P, 0] , which are created by explosive ionzation of a solid 
target by an incident laser. Fields comparable to stellar 
magnitudes (> 300 MG) have already been experimen- 
tally realized in the laboratory with the use of super- 
intense ultra-short laser pulses 3'j. Recent studies have 
also established the ultrashort duration picosec.) of 
these self-generated magnetic fields under femtosecond 
laser excitation 0, The temporal behaviour of the 
magnetic pulse contains crucial information about the 
dynamics of hot electron propagation inside a solid den- 
sity material / plasma. This information about the pen- 
etration and transport of hot electrons 0] inside an over- 
dense plasma and factors affecting it, is of utmost im- 
ortance to the fast ignitor approach of fusion research 



The pulsed high magnetic fields generated primarily 
near the critical density surface, result from a combina- 
tion of currents, viz. the direct hot electron current which 
is due to copious generation of fast electrons at the criti- 
cal layer via coUisionless absorption of incident laser en- 
ergy , and the return shielding current which is generated 
in response to the hot electron current [^. The propa- 
gation of hot electrons inside the solid density target de- 
pends on whether the target is a dielectric or a conductor. 
In case of a dielectric ( material with high classical resis- 
tivity ), the return shielding current is weak, resulting in 
a large charge imbalance which in turn exerts a strong 
retarding field on the fast electrons, resulting in electro- 
static inhibition. Electrostatic inhibition of hot electrons 
in dense matter has been seen in some experiments by 



measuring K a emission from layered targets. In initial 
experiments by Bond et. al. JlCll inhibition of fast elec- 
trons due to resistive electric field was observed in low 
density gold targets. In contrast to a dielectric, in a con- 
ductor ( material with extremely low classical resistivity 
), the return shielding current nearly balances the hot 
electron current. Here the stopping of hot electrons oc- 
curs because of turbulence induced anomalous resistivity, 
the chaotic magnetic turbulence effects being generated 
by t he p resence of large gradients in the plasma currents 
|l2l |l3j- Experimental evidence of turbulence induced 
anomalous resistivity is presented in a recently published 
experiment by Sandhu et. al. where they irradiate 
a solid aluminium target with a high intensity femtosec- 
ond laser. Since the mechanism of hot electron transport 
( along with return shielding current ) and their eventual 
stopping is different ( it is either electrostatically induced 
jo] or turbulence induced ) in a dielctric and a con- 
ductor, it is expected that the temporal evolution of the 
magnetic fields which these currents generate will also be 
different. Thus the complex dynamics of the hot electron 
propagation inside a solid density material ( or plasma ) 
is reflected in the time evolution of the magnetic pulse. 

Here we present experimental results on transport of 
hot electrons and their inhibition in both high and low re- 
sistivity media via measurements on temporal evolution 
of megagauss magnetic pulses. Recently, there has been 
a lot of experimental work on the transport of hot elec- 
trons but to our knowledge no experimental result 
exists which gives a clear distinction between electrostat- 
ically induced and turbulence induced inhibition of fast 
electrons. In our experiment, this has become possible 
through the use of temporal measurement of magnetic 
pulses which is intimately related to the propagation of 
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FIG. 1: Experimental setup, Target (T), Lens (L), Beam 
Splitter (BS), Photo Diode (PD), Half Wave Plate (HWP). 
Inset represents pump and probe paths in plasma. Expanded 
view of interaction region shows magnetic field geometry. 



hot electrons inside the material. We compare two en- 
tirely different media, one conducting (sohd Aluminum) 
and the other insulating (BK7 Glass). We use pump- 
probe polarimetry to measure the ellipticity produced 
in a probe laser pulse due to pump-generated magnetic 
fields. This measured ellipticity is then analyzed numer- 
ically to deduce the magnetic field as a function of time 
delay between pump and probe pulses. We find that 
the evolution of magnetic field in aluminium and glass 
to be drastically different. The time required for the 
magnetic field to decay in glass is found to be an or- 
der of magnitude smaller than in aluminium. This is a 
clear indication of electrostatic inhibition of fast electrons 
which results from absence of return shielding currents 
due to very low electrical conductivity of glass. We also 
present an one-dimensional model ( an improvement over 
our zero-dimensional model 4] ) for magnetic field evo- 
lution which gives an interpretation of our observations. 
The decay time in each case is related to the electrical 
conductivity ( electrostatically induced or turbulence in- 
duced ), thus providing us with estimates of conductivity 
of glass and aluminium under conditions of elevated tem- 
peratures and pressures. 

The experimental set-up is shown in figure 1. The 
p-polarized pump laser (lOOfs, 800nm) is focused at an 
angle of 45°, to an intensity of ~ 2 x lO^^Wcm^^ on a 
spot size of 20 microns. The long term contrast for the 
pump is 10^. The probe laser used in measurements is 
frequency doubled (400nm) in a thin BBO crystal, and is 
made ~ 10'^ — 10^ times weaker than the pump intensity 
and is focused at near normal incidence. In this configu- 



ration, the probe penetrates beyond the critical density 
(ric) for the pump laser (i.e. upto inf^^^), and thus 
samples high density region near and beyond the critical 
layer, where high magnetic fields are expected to exist. 
It should be noted that this high density region is not 
accessible using tangential polarimetric measurements of 
Faraday rotation. The input laser fluctuations are moni- 
tored with a photodiode (PDl). The reflected probe from 
the interaction region is split into two arms, one arm is 
used in PD2 to measure the reflectivity of the probe and 
the second arm is used in the quarter wave plate and 
polarizer combination ( kept in front of photodiode PD3 
) to measure the polarization parameters of the probe. 
The reflectivity measurement serves to establish the zero 
delay and is used to normalize other polarimetric mea- 
surements. All the measurements are made as a function 
of time delay between pump and probe. Diffuse scatter 
is also monitored with another PD and is observed to be 
insignificant and constant in our temporal range of inter- 
est. In our experiment, we have used two type of targets, 
first one an aluminium target i.e. a conductor and second 
one a glass target i.e. an insulator. The BK7 glass target 
was coated with about half a micron thick aluminium on 
the front surface so that laser hitting the target generates 
identical hot electron temperatures in both cases. To es- 
tablish that this is indeed so, we have measured the hard 
X-rays arising from hot electrons in the energy range from 
50 to 500 keV for both the cases, under similar laser con- 
ditions with a Nal(Tl) detector. The hot electron spectra 
are shown in figure 2(a). The hot electron temperature 
obtained from the fit (shown as solid line) is observed to 
be approximately same (~ 30fceT^) in both the cases (Al 
layer on Glass and Al metal alone). This indicates that 
we have identical hot electron source in both cases arising 
from interaction of the laser pulse with the front surface 
of the target. The magnetic field induced ellipticity in 
the probe pulse is measured using different combinations 
of quarter wave plate and polarizer in analyzer setup i.e. 
polarizer at 0, 45, 90° and A/4 plate at 45° and polarizer 
at 90° The measurement of stokes parameters in 

this method rules out the presence of any random de- 
polarization. Experimental curve showing ellipticity vs 
time delay is also shown in figure 2(b). The significant 
difference is observed in time scale when we compare the 
two cases, that of glass and aluminum. We now numeri- 
cally analyze the measured ellipticity in order to deduce 
the magnetic field present in the interaction region. 

An electromagnetic wave travelling through a magne- 
tized plasma with wave vector perpendicular to the qua- 
sistatic magnetic field can be analyzed in terms of two 
characteristic modes, the ordinary mode (0-wave) with 
electric field polarization vector parallel to the magnetic 
field and the extraordinary mode (X-wave) with electric 
field polarization vector perpendicular to the magnetic 
field l2l|. These two modes have different refractive in- 
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FIG. 2: The hard X-ray emmission spectra (a) for glass sam- 
ple coated with 0.6 micron Al (b) pure Al metal target; in en- 
ergy range 50-500 keV with laser intensity of 2 x 10^^ Wcm~^. 
Solid line is the exponential fit to deduce the hot electron tem- 
perature in each case, (c) Raw data for ellipticity in case of 
Glass coated Al (d) in case of solid Al metal. 



dices, which are 
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where /io represents refractive index for O-mode and 
represents refractive index for X-mode. Here uj, ujpe 
and LOce are respectively the laser frequency, the elec- 
tron plasma frequency, and the electron cyclotron fre- 
quency i.e. eB/mc. The difference in refractive indices 
results in accumulation of phase difference between O- 
and X- waves as they propagate through the plasma. 
This difference in phase is observed as induced ellip- 
ticity. If we represent the laser polarization in terms 
of stokes vector s = (si, 32,53) , si = cos{2x)cos{2'4>), 
S2 — cos{2x)sin(2ip) , S3 = sin{2x), where is ori- 
entation of polarization ellipse, then the ellipticity e is 
given by e = tan{x)- The evolution of stokes vec- 
tor inside the magnetized plasma is given by the differ- 
ence in refractive indices of the characteristic waves i.e. 
ds/dz — Q,{z) X s{z), where Q. = {uj/c){j2o — A'x)- The 
solution of this polarization evolution equation is writ- 
ten as s{z) — M.Sinput, where Mij is the transition ma- 



trix which obeys the e quat ion dMi/dz 
= {Mu,M2^,M^,) M 



VL X Mi, where 
We integrate this equation 



numerically inside the plasma by dividing it into small 
slabs, where within each slab, the plasma parameters are 
assumed to be constant. 

For the purpose of integration, the required plasma 
density profile is modelled by assuming that the plasma 
expands into vacuum at ion sound speed in a self sim- 
ilar fashion 0, [i^l, which yields an exponential den- 
sity profile of the form Ue = nsexp{— z / Cgt) , where n,, 
is the solid density and Cs — \/ (ZkT^jmi) is the ion 
sound speed. We have chosen the background electron 
temperature Tj, = lOOeF and the average ionization 
state Z = h [22|. At the vacuum end, which is cho- 
sen at a distance La from the critical layer, we put a 
cut off on the density ^n^nd which is chosen sufficiently 
low so that it makes practically no difference to the re- 
sults. The plasma slab to be integrated over is given by 
hs = Cst\ln(ncriticailnendf\- It IS found that the results 
are not sensitive to the chosen density profile; even the 
use of linear density profile yields results within accept- 
able error bars. In addition to the density profile, we 
have assumed uniform magnetic field inside the plasma 
column. Although, this assumption is not strictly cor- 
rect, it makes little difference to the results, and is a 
reasonable one as shown later. The solution Mij{z), and 
s{z) = M.Sinput can now be determined for any small 
plasma slab ( where the density is assumed to be con- 
stant ) for a given value of magnetic field B in definition 
of fl. Thus recursively solving in small steps and main- 
taining unitarity of M one can find the full transition 
matrix for whole plasma column, for any chosen density 
profile. From the transition matrix, the final stokes vec- 
tor can be computed which in turn gives the ellipticity 
of the reflected probe. 

In addition to the phase difference introduced by dif- 
ferences in refractive indices, there is an extra phase dif- 
ference which aries due to slightly different turning points 
for O and X-wave. The X-mode has a turning point i.e. 
right cut-off in the dispersion diagram, earlier than the 
O-mode, which reflects at usual critical density. This ex- 
tra path travelled by O-mode as compared to X-mode 
leads to an additional phase difference, which we have 
included in our calculation of net ellipticity. Thus the 
effect of magnetic field manifests itself in the form of el- 
lipticity of the reflected probe light by both of the above 
mentioned processes i.e. difference in refractive indices 
and difference in turning points. 

Using the above described procedure, for each time de- 
lay, we iteratively calculate the magnetic field B required 
to yield the experimentally observed ellipticity. The fi- 
nal plots of magnetic field for solid aluminium target and 
BK7 glass coated with aluminium layer are shown in fig. 
3(a) and fig. 3(b) respectively. The squares represent 
experimental points and the solid curves, which show 
a reasonable fit to the experimental data are obtained 
using a one-dimensional model developed below. It is 
observed that the peak magnetic field reached in both 
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FIG. 3: Magnetic field pulse profile for (a) Glass (b) Alu- 
minum with laser with intensity of 2 x 10^® Wcm~^. Solid 
line shows the fit for p-polarized case using a phenomenolog- 
ical model. 



cases are of the same order ( 32MG in case of BK7 glass 
and 22MG in case of solid aluminium ). The most strik- 
ing difference between the two cases is seen in the time 
durations of the self-generated magnetic pulses (note dif- 
ferent time axis of fig 4(a) and (b)). In case of BK7 glass 
coated with alumunium, the time duration of the mag- 
netic pulse ('^ Ipicosec) is an order of magnitude smaller 
than the time duration (~ lOpicosec) of the magnetic 
pulse for pure aluminium. This is a clear indication of 
high background resistivity of glass, which curtails the 
return currents and hence results in inhibition of hot elec- 
trons because of electrostatic fields. This leads to much 
faster resistivite decay of magnetic fields in case of glass 
sample as compared to solid aluminium sample, where 
the decay is mainly due to turbulence induced anoma- 
lous resistivity ,4]. 

We now model the temporal evolution of magnetic 
field using the following equation, which describes the 
mechanism of quasistatic magnetic field generation un- 
der EMHD approximation. 



as 

'dt 



Ana 



V^B + - (V X jhot 

cr V 



(3) 



The above equation is derived by taking the curl of the 
equation of motion of background plasma electrons that 
carry the return shielding currents, jp = (c/47r)(V x B) — 
jhot, where jp = aE, a being the conductivity of the 
background plasma. The first term describes the mag- 
netic field decay due to resistive damping of the plasma 
shielding currents and the second term describes the mag- 
netic field generation due to hot electron currents. As- 



suming the magnetic field to be in the azimuthal direc- 
tion, equation 10) can be approximated as 



'dt 



Brf, 



Siz,t) 



(4) 



where the diffusion term is approximated as B/t with 
T w (47rcr/c^)(Ar)2, and the source term is approximated 
as S{z,t) « —{c/(T{Ar))jhot{z,t). Here Ar is the laser 
spot radius ~ 10/i. Our main interest lies in understand- 
ing the evolution of the magnetic field after the laser 
is switched off. This is because the experimental results 
show that, although the laser is on only for ^ 100 f msec, 
the magnetic field peaks and remains on for picoseconds 
time duration. Taking B — Bias at t — Tiaser, the solu- 
tion of equation Q is given by 



B = Biase-^*-^'--^^^ + e~'/^ 



aAr 



Jhot 



where jhot = —s-rihVh, is the hot electron density and 
Vfi is the velocity of the hot electron fluid. To make an 
estimate of rih and Vh, we use the formalism given by 
Bell et. al. 0|. According to 0, the evolution of 
hot electron density {uh) is governed by the following 
nonlinear diffusion equation 



dt 



d ( aTh drih 



dz Ke^Uh dz 



(6) 



Here Th is the hot electron temperature. It can be shown 
that the above equation is valid even without the re- 
stricted assumption j total — jhot + jp ~ 0, used in ref. 

Since our interest lies in t > Tiaser, we use the solu- 
tion of equation © in this temporal regime, as given in 
ref. El 



rih 



2noZo 



z^ + L"^ 



with 
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2 labs'^laserS- 



6s 



(7) 

(8) 
(9) 
(10) 



where the absorbed intensity labs = flincident, f being 
the fraction absorbed. Here no is the density of hot elec- 
trons at z = 0, at time t — Tiaser and zq is the charac- 
teristic stopping length such that uoZq is the total num- 
ber of hot electrons produced at time t = Tiaser- The 
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constants no and zq have been derived by equating the 
absorbed laser energy ^%,bsTiaser" to the hot electron ki- 
netic energy (see ref. The above solution {71) is a 
self-similar solution of equation © in which the spatial 
shape remains the same but it expands in time with a 
scale length L{t). Using the above expression for Uh and 
L(t), we estimate jhot as 




Jhot 



-euhVh 



(a) 



I 



2nQZQ L 



dL 
dt 



(11) 



where Vh, the hot electron velocity is taken to be pro- 
portional to a being the proportionality constant. 

Substituting the expression for ^ in jhot and using it in 

2 

equation JSjl along with a « 4^(Ar)^ ^® 



B 



Ae 



z2 + z2(pyr-H 1)6/5 



dy (12) 



2^f^ andp 



We 



where w = ^^LJU^^ilI = — ■ ■■ ^^^^ „ = ^ — ^_ 

now use the above expression for B(t, z) at z = to 
model the magnetic field evolution as a function of time, 
for both BK7 glass coated with aluminium and solid alu- 
minium using r ( which is related to conductivity a ) and 
/ ( fraction of light absorbed ) as free parameters, z = 
is the point where hot electrons are generated; so it is 
actually the critical density point. The proportionality 
contstant a is taken to be unity. Bias , the magnetic field 
at time t = Tiaser is 15MG and 18MG respectively for 
aluminium and glass. The best fit curves ( solid lines ) 
are shown in fig. 3(a) and fig. 3(b) respectively. From 
the best fit, the relevant parameters ( / and t ) for alu- 
minium and glass are: 

Aluminum: f = 0.5; r = 3.5 ps 

Glass: f = 0.2; r = 0.2 ps 

Before presenting a discussion of our results obtained 
from one-dimensional modelling and numerical fit, we 
first justify the assumption we have made in our numer- 
ical analysis of ellipticity viz. the assumption of uniform 
magnetic field inside the plasma slab. Fig. 4 shows the 
typical magnetic field as a function of z and t obtained 
from our one-dimensional model. This plot is for alu- 
minium with / and r as calculated above. As expected, 
the magnetic field diffuses in z and decreases in ampli- 
tude as time increases. Although the magnetic field B 
varies with z, we shall see below that its variation is not 
significant over the range of scale lengths which are rele- 
vant for our experiments. In Fig. 5, using equation 1)12(1 
we have plotted the magnetic field as a function of posi- 
tion at a time corresponding to the peak of the magnetic 
field, for both aluminium and glass. In this plot, we have 
used the same value of / and r as obtained above. This 
plot clearly shows that the variation of magnetic field 
over the plasma slab is negligible. This justifies that the 
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FIG. 4: (a) The z and t dependence of the magnetic field 
profile inside the plasma based on one dimensional model for 
propagation of hot electrons as disscused in the text.(b) Spa- 
tial dependence at the time corresponding to the peak mag- 
netic field in both the cases. 



assumption of uniform magnetic field in our analysis of 
ellipticity is a reasonable one. 

We now discuss our results. From the results of the 
fit, it is clear that the time for magnetic field decay in 
aluminium is an order of magnitude larger than in glass. 
Calculation of conductivity using a = (c^/47r(Ar)^)r, 
gives GAi = 2.5 X 10^'^s"^ and Ogiass = 1-43 x lO^^s"^ As 
of now there are no measurements of conductivity at ele- 
vated temperatures for glass and very few in aluminium. 
We see that resistivity of aluminium as deduced from 
these measurements is an order of magnitude higher than 
the value of resistivity deduced from reflectivity measure- 
ments by Milchberg et. al. This, we believe, is due 
to turbulence induced anomalous resistivity which affects 
the return currents. This conclusion is also supported by 
simulation where similar effect of anamalous resistivity 
and stopping of hot electrons have been seen . In the 
case of glass, the situation is more complex as at low 
temperatures glass is almost non-conducting (cr ~ 0). 
However, at elevated temperatures, as a result of tar- 
get heating by coUisional effects and large electric fields 
exceeding breakdown threshold, conductivity becomes fi- 
nite. There exists observation of 8/i wide collimated ion- 
ization tracks extending upto 150 — 300/i and lasting for 
60 ps in glass with a 10^'^ Wcm~^ intensity laser Al- 
though we cannot estimate here the amount of target 
heating or ionization, we do get an estimate of resitivity, 
which we assume to be uniform in space and time. Our 
measurements show that neutralization of hot electron 
current is clearly not as effective in glass as it is in alu- 
minium, due to high background resistivity. This results 
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in inhibition of hot electron propagation in glass through 
generation large of electrostatic fields. We now make 
an estimate of the penetration depth (Lj) of hot elec- 
trons using equation (8) . Using the parameters Zq and a 
from the fit and Thot from observations, we get, for glass, 
i/ = 4.7 X 10~^cm in Ips and 2.1 x 10~*cm in lOps, and 
for aluminium Lj- = 1.6 x 10~^ cm in lOps. These num- 
bers indicate that hot electrons are penetrating an order 
of magnitude more distance in a conducting background. 
Similar results on electrostatic inhibition have also been 
measured in plastic targets using Ka emission from lay- 
ered targets. This technique has its limitations due to its 
reliance on x-ray measurements which is a secondary di- 
agnostic (with limitations of sensitivity, calibration, and 
reabsorption) and on numerical codes for interpreting X- 
ray data. 

In conclusion, the penetration of fast electrons into a 
target can result in very high electric field, which is not 
quenched by thermal return currents. This can lead to 
severe impediment in the flow of fast electron currents. 
We provide evidence of inhibition of fast electrons in an 
insulating media. Our results on magnetic pulse mea- 
surement and modelling in terms of hot electron currents 
as source, yield measurements of conductivity of conduct- 
ing and dielectric media under extreme conditions. These 
results are very important for laser fusion schemes, which 
rely on ignition of hot spot by energy deposition by fast 
electrons. 
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